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Abstract 

The effect of solar magnetic activity on the yearly mean average temperature is 
extracted from the historical record for much of North America. The level of solar 
activity is derived from the international sunspot number by the renormalized 
continuous wavelet transform using the Morlet basis to provide a running estimate 
of the power associated with the magnetic cycle. The solar activity gives the 
abscissa for a scatter plot of temperature for each station, from which the solar 
dependence and mean temperature are extracted. These parameters are then 
plotted against the latitude, longitude, and elevation for each station, revealing a 
dependence of their values on geophysical location. A mechanism to explain the 
latitudinal variation of the solar dependence is suggested. 



1 Introduction 

The amount of influence which the solar magnet i c cycle has on Earth ' s climate continues 



to be a mat t er of debate flLean and Rindl. 



2006 



2007: Lockwood and Frohlich 



9981: Wagner and Zo 



120071 lEcher et all . [2009; 



it a . 



2005 



Li et a 



Moore et all . 



20091: 



Sitnov 



2009). The intriguing coincidence of the Maunder Minimum (lEddyi . ll976uLuterbacher et all . 
200l[ ) , a period of exceptionally low solar activity as indicated by a deart h of sunspots 



observed fro m 1645 and 1715, and the coldest years of the Little Ice Age ( iGrovd . 12001 



Mannl . 120021 1 . a period o f exceptionally low temperatures as observed in Europe, North 
America, and elsewhere ( iKreutz et all . 119971 : iJohnson et all . 12001c iHolmgren et all . 1200 ll ) 
ending in the middle of the 19th century, suggests that the temperature on Earth 
might display a dependence on the level of solar magnetic activity. Elucidating that 
dependence is this article's primary topic of investigation. 
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Using the renormalized continuous wavelet transform to provid e a running esti- 



mate of the solar activity from the international sunspot number (jJohnsonl . [2009a 



2010al ). here we extend its comparison to the yearly mean average te mperature pro- 



vided by the United States Historical Climat ology Network, or USHCN (jOuinlan et all . 



19871 ; iKarl et all . Il990l ; lEasterling et all . 119961 ). for 1218 stations covering much of North 
America. The mean average temperature series as given does not come with an error 
estimate (though one may be derived from the underlying monthly measurements and 
is the subject of further investigation), thus one unit of variance is assumed through- 
out. From these data we extract the solar dependence and mean temperature for each 
station. 

Our basic strategy is to replace the abscissa of time in the temperature record for a 
station with the corresponding level of solar magnetic activity, thereby determining from 
the parameters of a linear regression that station's temperature's dependence on solar 
activity as well as the mean temperature without regard to the activity level. Those 
parameters then become the ordinate against the three abscissas of station coordinates. 
The distribution of these parameters displays a correlation with the geophysical location 
of the station, particularly for latitude and elevation but not so much for longitude, 
and another linear regression reveals their dependence. The latitudinal variation of the 
solar dependence might be explained by an increased polar influx of energetic particles 
from the solar wind during periods of greater solar magnetic activity. 



2 Solar magnetic activity 



To evaluate the level of solar activity from the historical sunspot record, we u se the 
renorm alized continuous wavelet transform. The algorithm is discussed in detail in I Johnson 
(l2009bl ). thus only some highlights are given here. Writing the Morlet wavelet at scale 
s > as the product of a real constant C s , a Gaussian window $ s = exp(— 1] 2 /2), and 
a normalized wave Q s = exp(iu>ir]) in terms of the parameter rj = (£' — t)/s, 



^ t {tf) = C s <S> Sjt (t')Q s , t (t') , 



(la) 
(lb) 



where u\ ~ 2tt is the central frequency of the mother wavelet at unity scale s = 1 and 
zero offset t = which here spans 2% + 1 = 13 time units At = 1, the normalization 
C s = (4/7TS 6 ) 1 / 4 produces a wavelet with norm \/2/s and a symmetric forward and 
inverse transform pair given a mean-subtracted signal y(t) —ty(t) — (y(t)) of duration 
N t , 



cwt(m) = £VC lt (Ov(0 



ICWT(t) 



Re 



E£^#)CWT( s ,i')A S 



s t' 



(2) 
(3) 



cf. Equations (6) and (9) by iFrick et all ( 119971 ). One compensates for the wavelet 



truncation caused by the finite signal duration by taking ip Btt ijj S: t\/2/s 2 \^ Sj t\ 2 which 
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preserves the wavelet norm. The peak normalized power spectral density PSD(s,t) = 
|v2CWT| 2 is twice the square of the amplitude of the CWT and displays in the instant 
wavelet power IWP 4 (s) = PSD(s,t) response peaks whose integrated area equals the 
sum of the squared amplitudes for an infinite signal with sinusoidal components of sta- 
tionary amplitude and period. The integrated instant power IIP(t) = J2 a PSD(s, t)As 
then gives a running estimate of the instant signal power, and the mean wavelet power 
MWP(s) = Nf 1 Y,t PSD(s, t) provides the net power spectrum. We have found satis- 
factory if not superior performance for data an alysis when neglecting the admissibil- 
ity condition ( Frick et al . 1997 ; Johnson! . 2010a ) related to the DC response (nonzero 
mean) of the analyzing wavelet J2t' V's.tC^O 7^ 0> which remains negligibly small until the 
wavelet truncation becomes significant. 

Since the invention of the telescope, solar astronomers have been able to keep an 
instrumental record of the magnetic flux erupting from the surface of the sun, which 
appears as d arkened spots. First for mulated by Wolf, the international sunspot num- 
ber Rj by thelSIDC-teaml (11700-20081 ) continues the Zurich number up until the present 
day (IHossfieldl. 120021) and is know n to display a tighter correlation with sunspot area and 
radio flux (IHathaway et all . l2002l ) than the sunspot group number R g by lHoyt and Schatten 
(119981 ). The CWT analysis of the mean-subtracted yearly smoothed sunspot number 
Ri is presented in Figure [lj where the PSD has indicated in white the cone-of-influence, 
where wavelet compensation (which begins at the more restrictive cone-of-admissibility) 
becomes significant, as well as the scale of the signal duration at 309 years, beyond which 
lies the extremely low frequency ELF region. Prominent peaks in the MWP are indi- 
cated by the dotted lines. Its IIP gives the spot number activity SNA and outlines the 
magnitude of Ri occuring at the maxima of the Schwabe cycle roughly every 11 years. 

Unde r lying the sunspot number cycle is the solar magnetic activity discovered 
by iHald (119081) w i th a period of about 22 years and extending over the entire solar 
surface (iBabcockl . Il96ll ). The magnetic field of the spots reverses polarity over one 
Hale cycle, thus the sunspot number represents in some sense a re ctified version of the 
solar magne tic activity signal. In accord with the investigation bv lBuck and Macaulay 
(1992, 119931 ) and after a thorough evaluation of the arguments by iBracewelll (11953 . 



1983), we apply alternating signs to each of the Schwabe cycles in Ri to form the derec- 



tified solar magnetic signal ±Ri. Its CWT analysis, shown in Figure^ displays a much 
simpler spectral content dominated by the Hale cycle and its third harmonic at about 
7 years. Its IIP gives the solar magnetic activity SMA and is approximately double 
the SNA, with either denoted b y SA. Both the SNA and the SMA compare well with 
the power estimate fo und by iLd ( 20041) de termined from just the Schwabe cycle. While 
aware of the work by lUsoskin et all (120091 ) . we maintain the traditional numbering of a 
tenth Schwabe cycle beginning about 1800. 



3 USHCN temperature record 

Unfortunately, the North American historical climate record does not begin shortly after 
the invention of the thermometer but gradually picks up contributing stations through- 
out the 19th and early 20th centuries. For this analysis, we utilize the USHCN yearly 
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mean average temperature data ( Quinlan et al . 1987 : Karl et al . 1990l : Easterling et al . 
19961 ) as our climate signal. The USHCN stations cover much of North America, as 
depicted in Figure El spanning approximately 30 degrees of latitude, 60 degrees of longi- 
tude, and 3 kilometers of elevation. The density and extent of the station data reflects 
the duration of Western occupation as well as the accessibility of the location. The 
values carry a flag indicating where missing data has been replaced with an estimate 
from surrounding data, and we will consider both the full data set as well as a restricted 
set consisting only of values flagged "good" . 

In performing our linear fits, we adopt a prior which is uniform on the angle of the 
slope rather than its magnitude, summarizing our methodology for those u nfamiliar 



with Bayes i an da ta analysis expressed in terms of conditional probabilities (IDurrett 
1994 ISivial . Il996l ). Using notation p(M\ I D) = prob(M|D,7) reading "probability of 
M given D and information 7", we further abbreviate p(M\iD) to p(M\D) = pj$ and 
p(M\I) = p M when the background information I is unchanging. For model M with 
parameters m = {rrij} and data D = {D^}, Bayes' Theorem allows one to write 
Pd = PmP M /p D reading "the posterior for M\D equals the product of the likelihood 
for D\M and the prior for M divided by the evidence for D" , where the denominator 
p D affecting neither parameter estimation nor model selection is often omitted. The 
best estimate for the model parameters rap is given by the maximum of the posterior 
oc PmP™ ', whose logarithm (base e) is written as L P = L L + L m — L D , where the 
log evidence L D = is a constant. 

For a uniform prior L m = =ff m , the maximum of the posterior is at the maximum of 

with Gaussian noise 



rriL, and for independent data pfy 



the likelihood mp 

PZ = UP% = Il(2^)- 1/2 exp(-^./2) , (4) 

k k 

where = (Mj, — D^j/a^ is the weighted residual for the fcth datum, so that mi 
minimizes the least-squares residual x 2 = J2k R\- F° r a linear model Mk = yk = 
mxk + c giving the ordinate over an abscissa x, the solution {m, c}i is found by setting 
-VL L = Vx 2 / 2 = °> where 



1. 



-yx 2 



5*20 


•S'lo 




m 




Su 


5*10 


Soo 




c 




5*01 



(5) 



in terms of the weighted sums Sab = J2k x t^ > k/ a ky an< i ^ one takes Xk —> x' k = Xk — (x), 
the parameter cl returns the best estimate for the mean (D) and the covariance matrix 
will be diagonal so that the error correlation r a = a 2 mc / ' (cr^^a^) 1 ^ 2 is zero, where the 
mean is computed with respect to weight J2k°~k 2 - Linearity implies that the shape of 
the posterior is a Gaussian given solely by the second derivative, 



VVL, 



1 9 

2 VV * 



5*20 


•S'lo 




2 2 


-1 


Sio 






2 2 
® mc ® cc 





(6) 



whose width determines the covariance matrix. With unit variance on the data, the 
deviation of the parameter s depends only upon the valu es of the abscissa, which we 
assume are known exactly (IPress et all . 11992c ISivial . Il996[ ) . 
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A prior which is uniform on the magnitude of the slope m will not be uniform on 
the angle of the slope 9. To show no preference for the slope of a line, one should take 
a prior which is uniform on its angle, and given unit bearing axes x and y, that angle 
should be determined on the scale invariant axes X and Y to be defined. First, one 
takes X = x'/A x for Ax = max(|a4|) to scale the abscissa to [—1,1], and then the 
ordinate Y = y/Ay for Ay = max.(\D^ — (D)\) is not centered but is scaled as if it 
were. On these axes, the line has slope m' = mAx/Ay with angle 9 = arctanm', and 
it is on 9 £ (— 7r/2,7r/2) that we apply the uniform prior p e = 1/ir to encompass slopes 
within the range m ~ mf G (— oo, oo). In terms of m', one writes p e = p m |dm'/d#| to 
find p m = [ir(l + m' 2 )] -1 such that — L m > = log(l + m' 2 ) + # n , and then in terms of m 
the prior is Lorentzian, 



V" 



A 



-7T 



X 



(7) 



with normalization / p m dm' = f p m dm = 1. For comparison, a prior uniform on m! 
will appear in 9 as p e oc 1 + (tan#) 2 , which clearly displays a preference for a slope of 
extreme magnitude. 

The prior nonuniform on m ~ m' introduces a nonlinearity to the posterior, which 
becomes Lp = Ll + L m + L c — Lp = Ll + L m + # c ,r>, as p c = 1/A C over a range of 
A c . To find m'p, one must solve VLp = 0, where 



VL, 



-Vx + Vlog(l +rr£ 



'2 \ 



10 



10 
ql 

°00 



m 
„/ 



ql 
ql 



01 



+ 



2m'/ (1 + m' 2 ) 




lb) 



where S' ab are the weighted sums and {mf, c'} are the model parameters on the scale 
invariant axes X and Y. The equation for d remains linear, leaving a cubic equation 
to solve for vn! . The appropriate real root is found between m' L and 0, and the solution 
is then scaled {m',c'}p — > {m,c}p, as is the covariance matrix — [WLp] -1 , which 
differs from -[VVL L ] _1 only by S' 20 -)> S' 20 + 2(1 - m /2 )/(l + m' 2 ) 2 . The nonlinear 
term introduces higher order corrections to the shape of the posterior, which remains 
very close to a Gaussian unless dominated by the prior. 

The net effect of the Lorentzian prior is to prevent one from overestimating the 
magnitude of the slope best supported by the evidence. The estimate for the mean 
remains unaffected. In the limit of poor data |VL m | ^ |VL^|, the best estimate for 
the slope mp is driven from tul towards zero. We can see the effect of the prior on a 
fit to three points with increasing variance in Figure HI where the maximum likelihood 
estimate is a dashed line. One should consider how often one sees maximum likelihood 
used for cases such as (c) and (d), where the latter exhibits insufficient evidence for a 
significant slope. 

After so much digression, let us look at some data. For each of the 1218 stations 
indexed by j we take its temperature record and replace the abscissa of time with the 
value for the solar activity of that year (normalized by 10 5 ). In Figure [5] we display the 
linear regression for two locations chosen at random using the SMA and the full data 
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set, and in Table [T] we give the corresponding parameters for all four combinations of 
SA and data set. The mean station temperature T , indicated by the heavy dot in the 
figure, is also the value of the line 

T^(SA) = r sA xSA' + T J (9) 

at the mean of SA and indicates the temperature expected at that location without 
regard to solar magnetic activity. The solar dependence T gives the correlation of 
the mean average temperature with the solar magnetic activity and for any particular 
station may be in either direction, as seen in Figure 

There remains to answer the question of how much evidence there is for a nonva- 
nishing solar dependence T. One cannot simply ask how good is the fit, as the only 
logical answer is, compared to what? Bayesian data analysis frames the question in 
terms of model selection using the ratio of posteriors P = p^/Pd = PbP B /VaV A — 
(p B / p A ) (Pb / Pa) > where here model A has the single parameter T and model B has 
two parameters To and T. With no prior preference for either model p B jp A = 1, that 
ratio becomes 

_ II p r B p? 0iT:B dT dT 
V D A IP? 0>A dT ' 

where a common factor of p T ° = 1/Ar has been cancelled, whose logarithm (base 10) 
we write Piq. For the remaining prior factor we restrict the range of the angle to 7r/4 
so that = 4p m (m — > T) by arguing that any child could draw a line through a given 
set of dots with a slope correct to within a slice of it, and in Table [2] we give heuristic 
descriptions of several prior ranges which could be assigned. The question essentially 
boils down to determining how much to penalize model B for having an extra parameter, 
and we feel that half a quadrant is sufficient, as without the prior factor model B would 
always be preferred. In the numerator we take the quadratic approximation (which 
is exact for model A), p B ~ 27rarrO'T T exp(— Xt q ,t I^)Pb evaluated at the posterior 
maximum {T , T}p, after verifying its accuracy numerically. In Figure Owe display Pio 
for all stations using the SMA and the full data set, and in Table [3] we give its mean 
and median for both SA and data sets. We admit a preference for the SMA on the basis 
of the polarity associated with the Hale cycle, and we state subjectively that there is 
sufficient evidence to ascribe physical significance to the solar dependence T. 



4 Results 

Collecting our parameters into two sets {Tq} and {T J }, we proceed by plotting the 
ordered pairs (T j , X{) and (T j ,X{) for i E {1, 2, 3} and X j = (LAT J , LONG- 7 , ELEV- 7 ) 
and perform another linear regression. With so many points contributing, the prior 
here makes not much difference, but we maintain its use in our evaluation. Beginning 
with the mean station temperature To, we plot its dependence against station location 
X J in Figure [7] for the SMA and full data set, and in Table H] are the results of the 
linear regression for all the combinations. The column labelled J (x 2 ) gives the rms 
weighted residual, ie the expected distance between the fit and a data point. We see 
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that T displays the tightest correlation with latitude, as confirmed by the magnitude of 
Pearson's r = rp. Its dependence on latitude and elevation is not unexpected, and its 
dependence on longitude is minor, most likely the result of correlation among the 
representing the surface of the Earth. The units (hm) of elevation are chosen so that 
its range is approximately equal to that of latitude, thus the relative dependence given 
by their slopes may be compared directly We see that the decrease in T attributed to 
1 degree of increasing latitude equates roughly to that attributed to 2 hm of increasing 
elevation and that the magnitudes of their rp are in about the same ratio. The choice 
of solar activity signal does not affect T , and the choice of full or good temperature 
data has little impact on the results. 

The solar dependence T also displays a dependence on location, which again is 
strongest for latitude. Figure M displays the SMA and full data set results, and Table 
gives the entire set of results. The station averaged T is consistently positive, and we 
note a paucity of negative values at the lowest latitudes surveyed. Dots for stations 
with little evidence for a slope are shaded darker and confined to a band around T = 0. 
The solar magnetic activity SMA results display a tighter correlation with location 
than those of the spot number activity SNA, and its units were chosen to make the 
visible range of T roughly the same as that of To so that their relative slopes may also 
be compared. In these units, the dependence of T on latitude is about a third the 
magnitude and opposite in direction to that for T , and its r P has been reduced also by 
about a third. Its dependence on elevation is a third of that for latitude, and its variation 
with longitude is negligible. The tighter correlation and smaller errors and residual 
indicate that the solar magnetic activity derived from the derectified sunspot number 
provides a more physical signal than the activity derived from the usual (rectified) 
sunspot number. The choice of temperature data set does impact the results, with the 
full data set indicating a greater solar dependence. 



5 Discussion 



Comparison of our work to others' is made difficult because so few investigators con- 
sider the derectified sunspot number in their analysis of solar magnetic activity and 
most employ a cross coherence analysis looking for cycl es common to both the his- 
torical sun s pot r e cord and various indicators of c limat e (IBarnston and Livezeyl. 11989 



20051 : iBarlvaeva et all . 120091 : ISitnov 



Zhao et all . 12004 iPiscaronoft et all . 12004 : IWeber . 
20091 ) . Our methodology here is quite different, eliciting a relationship between the 
net solar magnetic activity derived from the instant power carried by the derectified 
sunspot signal and the temperature recorded across much of North America for almost 
a century and a half. Wh at that data lacks in du ration compared to the Central Eng- 
land Temperature record (jJohnsonl . l2009al . |2010aJ) it more than makes up for in breadth 
of ge ophysical location. We note t hat our results are in accord with othe r inves tiga- 



tors ( IFriis-Christensen and Lassenl . Il99ll ; iLean and Rindl . Il998l ; lEcher et all . 120091 ) and 
contra st with those who h ave found negligible evidence of sola r forcing on Ea rth's cli- 



mate flMoore et all . 12009 . 120071 : lLockwood and Frohlichl . 120071 : iLi et all . 120091 ). Most 



authors, when speaking of the influence of solar forcing on Earth's climate, are actually 
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discussing the e ffect of changing some model parameters in a simulation of meteoro- 
logical systems (jShindell et all . 120031 ; IWagner and Zorital . 120051 ). whereas here we are 
discussing a correlation found between independent instrumental observations from the 
historical record. 

While several mechanisms have been prop osed to exp l ain the possible coupling be- 
tween solar activ i ty an d terrestrial climate ( jSvensmarkl . Il998l ; iLean and Rindl . [1998; 



Hameed and Led . 120051 ). notably variation in irradiance, modulation of cosmic ray in- 



flux, and changes in global circulation, we would like to consider an alternative explana- 
tion. From the dependence of T on latitude, we see that the degree of solar dependence 
increases with distance from the equator. Thus, any proposed model for coupling should 
account for the greater dependence on solar magnetic activity towar ds the polar regions . 
One effect of increased solar activity is an increased solar wind (IWebb and Howard! . 
19941 ). whose temperature 10 5 K~ lOeV is sufficient for ionization to the plasma state. 



Those charged particles follow the lines of Earth's geomagnetic field emanating from 
the poles to produce the auroral illuminations and heating of the upper atmosphere. 
By the well known conservation law, that energy must be deposited somewhere, and 
we propose that during times of increased solar activity, an increased influx of energy 
bearing solar wind particles produces a preferential heating towards the polar regions. 
Whether that prediction will be supported by more detailed evaluations remains to be 
investigated. The first step in such an evaluation would be to verify the relation between 
the energy flux from the s olar wind and the solar magnetic activity, such as one can do 
for the 10.7cm radio flux ( Uohnsonl . l2010bl ). to put the historical sunspot record onto 
an axis bearing units relevant to geophysical analysis. Then, with an estimate of the 
historical thermal flux variation, one could look for correlations with various climate 
indicators. 



6 Conclusions 

A running estimate of the net solar magnetic activity is derived from the integrated 
instant power of the historical sunspot record, derectified to account for the opposite 
polarities of adjacent Schwabe cycles, using the renormalized continuous wavelet trans- 
form. The net solar magnetic activity provides the abscissa for a linear regression using 
a Lorentzian prior of temperatures observed at stations covering much of North Amer- 
ica for a duration exceeding one century. A relationship is found between the level of 
solar activity and the observed temperature characterized by the slope of the linear 
regression which itself depends upon geophysical location. The solar dependence is 
found to increase with latitude, suggesting a coupling mechanism to the polar regions 
of thermal flux from the solar wind. 
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Table 1: Linear regression results for two random stations 



ID# 


data 


SA 


T T (Ttt °"t t -Pio 


253630 




SMA 


7.2 47.7 2.01 0.0937 1.88 


SNA 


13.6 47.7 4.94 0.0937 0.826 


O 
O 
bC 


SMA 


5.32 47.7 2.12 0.0995 0.477 


SNA 


9.52 47.7 5.22 0.0995 -0.0718 


313017 




SMA 


-9.46 61.5 2.00 0.0937 4.11 


SNA 


-20.6 61.5 4.93 0.0937 3.12 


o 
o 

bO 


SMA 


-9.82 61.6 2.17 0.105 3.72 


SNA 


-20.8 61.6 5.31 0.105 2.72 



Table 2: Heuristic description of various prior ranges 



prior range 



the scattered dots are 



71 
7T/2 

tt/3 
tt/4 



there 

going up or going down 
going up, holding steady, or going down 
sharply up, slightly up, slightly down, or sharply down 



Table 3: Mean and median P 10 



data 


SA 


mean P 10 


median P 10 


§ 


SMA 


2.9172 1.0178 




SNA 


2.073 


0.51658 


o 


SMA 


1.6315 


0.19813 


o 

bC 


SNA 


1.0944 


-0.02401 
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Table 4: Linear regression results for T 





data 


m 


c 


^rom 


<y cc 


v^ 2 ) 


r P 




full 


-1.48 


52.4 


0.000539 


0.00269 


4.01 


-0.879 




good 


-1.48 


52.3 


0.000636 


0.00314 


3.86 


-0.884 


o 


full 


0.0562 


52.4 


0.000183 


0.00269 


8.37 


0.0982 


o 


good 


0.0526 


52.3 


0.000216 


0.00314 


8.22 


0.0927 


H 


full 


-0.637 


52.4 


0.000484 


0.00269 


7.63 


-0.42 


i-l 
H 


good 


-0.632 


52.3 


0.000588 


0.00314 


7.53 


-0.41 



Table 5: Linear regression results for T 





data 


SA 


m 


c 


@rnm 




\Ax 2 > 


r P 






SMA 


0.441 


4.52 


0.0116 


0.0576 


6.75 


0.31 


LAT 




SNA 


0.786 


9.3 


0.0284 


0.142 


14.9 


0.255 


o 


SMA 


0.442 


3.28 


0.0148 


0.0723 


7.35 


0.282 




o 


SNA 


0.678 


6.4 


0.035 


0.172 


15.8 


0.207 


•NG 




SMA 


-0.03 


4.52 


0.00392 


0.0576 


7.08 


-0.0621 


SNA 


0.0235 


9.3 


0.00965 


0.142 


15.4 


0.0224 


O 


o 


SMA 


-0.00781 


3.28 


0.00504 


0.0723 


7.66 


-0.0146 




o 


SNA 


0.1 


6.4 


0.0119 


0.172 


16.1 


0.0897 


> 
H 




SMA 


0.156 


4.52 


0.0104 


0.0576 


7.04 


0.122 


SNA 


0.241 


9.3 


0.0256 


0.142 


15.3 


0.0869 


i-J 

H 


o 


SMA 


0.142 


3.28 


0.0141 


0.0723 


7.63 


0.0952 




o 


SNA 


0.136 


6.4 


0.0329 


0.172 


16.1 


0.0443 
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Figure 1: CWT analysis of the international sunspot number whose IIP gives the spot 
number activity SNA 
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Figure 2: CWT analysis of the derectified international sunspot number whose IIP 
gives the solar magnetic activity SMA 
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Figure 3: Geophysical location of the USHCN stations 




Figure 4: Effect of the Lorentzian prior (solid) compared to the maximum likelihood 
estimate (dashed) as the variance increases 
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Figure 5: Linear regression of mean average temperature T against solar magnetic 
activity SMA using the full data set for two stations chosen at random 
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Figure 6: Logarithm (base 10) of the relative probability of the linear model for all 
stations using the solar magnetic activity SMA and full data set with its mean and 
median values marked 
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Figure 7: Linear regression for the mean station temperature To against geophysical 
location for the full data set 
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Figure 8: Linear regression for the mean solar dependence T against geophysical loca- 
tion for the solar magnetic activity SMA and full data set 
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